A qubit strongly-coupled to a resonant cavity: asymmetry of the spontaneous 
emission spectrum beyond the rotating wave approximation 
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We investigate the spontaneous emission spectrum of a qubit in a lossy resonant cavity. We 
use neither the rotating-wave approximation nor the Markov approximation. The qubit-cavity 
coupling strength is varied from weak, to strong, even to lower bound of the ultra-strong. For 
the weak-coupling case, the spontaneous emission spectrum of the qubit is a single peak, with its 
location depending on the spectral density of the qubit environment. Increasing the qubit-cavity 
coupling increases the asymmetry (the positions about the qubit energy spacing and heights of 
the two peaks) of the two spontaneous emission peaks (which are related to the vacuum Rabi 
splitting) more. Explicitly, for a qubit in a low- frequency intrinsic bath, the height asymmetry of 
the splitting peaks becomes larger, when the qubit-cavity coupling strength is increased. However, 
for a qubit in an Ohmic bath, the height asymmetry of the spectral peaks is inverted from the same 
case of the low-frequency bath, when the qubit is strongly coupled to the cavity. Increasing the 
qubit-cavity coupling to the lower bound of the ultra-strong regime, the height asymmetry of the 
left and right peak heights are inverted, which is consistent with the same case of low-frequency 
bath, only relatively weak. Therefore, our results explicitly show how the height asymmetry in the 
spontaneous emission spectrum peaks depends not only on the qubit-cavity coupling, but also on 
the type of intrinsic noise experienced by the qubit. 

Keywords: spontaneous emission spectrum, vacuum Rabi splitting, qubit strongly-coupled to a 
cavity 
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I. INTRODUCTION 

Strong and ultra-strong qubit-cavity interactions have 
been achieved in both cavity QED and circuit QED sys- 
tems (see, e.g., [H-01)- This opens up many new possible 
applications. For example, one could use the cavity as 
a quantum bus to couple widely-separated qubits in a 
quantum computer 0, [y] , as a quantum memory to store 
quantum information, or as a generator and detector of 
single microwave photons for quantum communications 
1- 

As a demonstration of strong interaction in cavity 
QED and circuit QED systems, the vacuum Rabi split- 
ting has been an exciting subfield of optics and solid- 
state physics (see, e.g., |7H10|). after its observation in 
atomic systems [ll|. In 2004, two groups [IH, re- 
ported the experimental realization of vacuum Rabi split- 
ting in semiconductor systems: a single quantum dot in 
a spacer of photonic crystal nanocavity and in semicon- 
ductor microcavity, respectively. In the same year, the 
experiment [Tij showed that the vacuum Rabi splitting 
can also been obtained in a superconducting two-level 
system, playing the role of an artificial atom, coupled to 



an on-chip cavity consisting of a superconducting trans- 
mission line resonator. When the qubit was resonantly 
coupled to the cavity mode, it was observed 14| that 
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two well-resolved spectral lines were separated by a vac- 
uum Rabi frequency i^R, a bi ~ 2g. Except for the asym- 
metry in the height of the two split energy-peaks (e.g., 
Ref. [3), the data is in agreement with the transmis- 
sion spectrum numerically calculated using the rotating 
wave approximation (RWA). When considering the vac- 
uum Rabi splitting behavior for strong qubit-cavity cou- 
pling, the anti-rotating terms should be taken into ac- 
count, and this might explain the observed asymmetric 
spontaneous emission (SE) spectrum. 



A. Antirotating terms are important for strong 
coupling QED 

It is obvious that the anti-rotating terms are not im- 
portant when the coupling between a qubit and the cavity 
field is sufficiently weak, and when the energy spacing of 
the qubit is resonant with the central frequency of the 
cavity. However, in the ultra-strong coupling regime, the 
anti-rotating terms of the intrinsic bath of the qubit cou- 
pling play an important role fT5T4l8l |. As a consequence 
of the anti-rotating terms in the Hamiltonian of cavity 
QED and circuit QED, even the ground state of the sys- 
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tern contains a finite number of virtual photons. The- 
oretical research [l9|, [2(| reveal that these virtual pho- 
tons can be released by a non-adiabatic manipulation, 
where the Rabi frequency g(t) is modulated in time at 
frequencies comparable or higher than the qubit transi- 
tion frequency. This phenomenon, called "emission of the 
quantum vacuum radiation" , would be completely absent 
if these anti-rotating terms are neglected. The energy 
shift of the qubit in its intrinsic bath has been studied 
in [2l[ using the full description, (i.e., non-Markov and 
without RWA) and found that the deviations from the 
previous approximation result already amount to ~ 5% 
for g/A ~ 0.1. 

In doped semiconductor quantum wells embedded in 
a microcavity (e.g., |22j]). considering the anti-rotating 
coupling of the intracavity photonic mode and the elec- 
tronic polarization mode, but using RWA in the coupling 
to their respective environments, it was found [22| that 
for a coherent photonic input, signatures of the ultra- 
strong coupling have been identified in the asymmetric 
and peculiar anticrossing of the polaritonic eigenmodes. 
From the descriptions given above, it can be seen that, 
as g/A increases to the ultra-strong coupling case, the 
anti-rotating terms that are otherwise negligible become 
more relevant and will lead to a profound modification 
in the nature of the quantum state of the qubit system. 



B. The asymmetry of the two splitting Rabi peaks 
can be explained beyond the RWA approximation 

In this paper, we study the SE spectrum of a qubit 
in a cavity. Our calculations include two kinds of anti- 
rotating terms: one from the intrinsic qubit environment 
and the other one from the cavity environment. This 
method is a powerful tool to investigate various kinds of 
qubit-environment interaction with anti-rotating terms 
and without using the Markov approximation. Because 
in this method, the qubit-environment coupling terms 
higher than the two-order are droped, it constricts that 
this method is unavailable in the case when the qubit- 
environment coupling str eng th is larger than the energy 
spacing of the qubit 0, [24[. Comparing the cases of a 
qubit in an Ohmic bath with the case of a qubit in a low- 
frequency bath, we find that for the case of a qubit in a 
low-frequency bath, as the qubit-cavity coupling strength 
increases, the height asymmetry of two splitting peaks is 
enhanced. However, for the case of a qubit in an Ohmic 
bath, the height asymmetry of the spectral peaks are 
inverted from the same case of the low-frequency bath, 
when the qubit is strongly coupled to the cavity. In- 
creasing the qubit-cavity coupling to ultra-strong regime, 
the height asymmetry of the left and right peaks are in- 
verted, which is consistent with the same case of low- 
frequency bath, only relatively weak. Since experiments 
reported that a superconducting qubit intrinsic bath is 
mainly due to low-frequency noise, our results are con- 
sistent with experimental data using a superconducting 



qubit in Ref. QJ]. 

We also investigate the dependence of the SE spec- 
trum on the strength of the qubit-cavity coupling and 
the quality factor Q of the cavity in either an Ohmic or 
in a low-frequency intrinsic qubit bath. Furthermore, we 
distinguish the contributions to the asymmetry from each 
bath: the intrinsic qubit bath and the cavity bath, and 
clarify the reason for the different kinds of peaks asym- 
metry. All of these results directly indicate that in the 
strong coupling regime, the SE spectrum is deeply influ- 
enced by the anti-rotating terms and the type of intrinsic 
noise experienced by the qubit. 

II. BEYOND THE ROTATING WAVE 
APPROXIMATION 

By using a cavity to confine the electromagnetic field, 
the strength of the qubit-cavity interaction can be in- 
creased by several orders of magnitude to the regime of 
strong or even ultra-strong coupling [25| . The strong- 
coupling regime for cavity quantum electrodynamics has 
been reached for natural atoms in optical cavities, su- 
perconducting qubits in circuit resonators (i.e., on-chip 
cavities) , and quantum dots in photonic-crystal nanocav- 
ities. Recently, the ultra-strong coupling regime has be 
achieved for a superconducting qubit in an on-chip cav- 
ity [2|. Although the coupling of the qubit to the cavity 
is much stronger than the coupling of the qubit to its 
intrinsic environment, the parameters in Ref. [14j show 
that both the decay rate of the cavity photon (k/2it rs 0.8 
MHz) and the qubit decoherence rate ('j/2tt « 0.7 MHz) 
are comparable. Therefore, we model the environment 
of the qubit in a cavity using two bosonic baths: one, 
called the intrinsic bath of the qubit and represented by 
operators b k and bl, is related to the relaxation of the 
qubit induced by its intrinsic environment; and the other, 
denoted as cavity bath of qubit and represented by the 
operators a\ and a k , involves the relaxation of the qubit 
caused by photons in the cavity. Figure 1 schematically 
shows the model considered here. For the intrinsic qubit 
bath, a broad frequency spectrum (e.g., either an Ohmic 
or a low-frequency spectrum) can be used to characterize 
it. For the cavity bath, because of the cavity leakage, it 
can be described by a Lorentzian spectrum with a central 
frequency, i.e., a single-mode cavity with its frequency 
broadened by the cavity leakage. The Hamiltonian can 
be written (throughout this paper, we choose h — 1) [22[ 
as 

H = -Aa z +^u kt iala k + ^g kj i(al + a k )o- x (l) 

k k 

+ y^a; fc|2 b[bfc + ^gk,2( b l + b k)o- x , 

k k 

where a x = a + + <r_ and i — 1,2, denote the intrinsic 
and cavity baths of the qubit, respectively. 

The baths experienced by the qubit can be character- 
ized by a spectral density Ji{uj) — J2k 9k — To 
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deal with the anti-rotating terms in Eq. ([I]), we will per- 
form a unitary transformation on the Hamiltonian. This 
unitary transformation is applied to the Hamiltonian H 
as follows: 



H 1 = exp(S)Hexp(-S), 
with S = Si + S2 , where Si (i — 1, 2) is given by 

S 2 = V* £ fc 2 (6+ - b k ) cr x . 

L — < Wh o. V / 



(2) 

(3) 
(4) 



In Eq. ([3j EJ, the parameter = uj k j/(uj k j +Vj A) is a 
fc-dependent variable. Up to order gr^, the transformed 
Hamiltonian ifW can be written as 

ij' pa I,y Acr 2 + ^2 Uk,ia\a k + ^2 u kt2 b\b k 

k k 

+ E 9k, 1 + a k cr + 



^9k,2 (b k <T- + b k <r + ^j , 



(5) 



where 



9k, 



2r/ 4 A 



+ r\i A 



v = mm, 



9k. 



r]i = exp 



E 



-7 



Sfc; 



(0) 



(7) 



(8) 




FIG. 1: (Color online) Schematic diagram of a two-level sys- 
tem or qubit with dissipation rate 7, which is coupled to a 
cavity with loss rate k by the qubit-cavity coupling strength 
9- 



and r\ = r\\r\i- Now the transformed Hamiltonian ([5]) has 
the same form as the Hamiltonian under the RWA, but its 
parameters have been renormalized to include the effects 
of the anti-rotating terms related with the intrinsic and 
cavity baths of the qubit. From the transformed Hamilto- 
nian iJ', one can see that, based on energy conservation, 
the ground state of the transformed Hamiltonian H' is 



\9 



(9) 



(using a z |4_) = — |4-)) and the corresponding ground-state 
energy is — 77 A/2. Therefore, the ground state of the 
original Hamiltonian H is given by 



|. 9 )= exp (-S) |s') 



(10) 



which is the dressed state of the qubit and the environ- 
ment due to the anti-rotating terms [l?J, [2?J] . 

A qubit can experience different types of intrinsic 
baths. The most commonly-used bath is the photon 
or phonon bath, which can be described by an Ohmic 
spectrum. However, for many solid-state qubits (e.g., 
superconducting qubits), the dominant dissipation can 
be due to two-level fluctuators, which behave like a low- 
frequency bath 26]. Here we consider either an Ohmic 
or a low-frequency intrinsic bath. The Ohmic bath with 
Drude cutoff is given by 



jp hm H 



2QOhmW 



1 + (cj/w hm) Z 



(11) 



where wohm is the high-frequency cutoff and aohm 
a dimensionless parameter characterizing the coupling 
strength between the qubit and its intrinsic bath. Here 
the low-frequency bath is written as 



j} ow H = 



2ai ow uj 



(uVA) 2 + (cj 1ow /A) 2: 



(12) 



where cji ow is a characteristic frequency lower than the 
qubit energy spacing A of the qubit, and ai ow is a di- 
mensionless coupling strength between the qubit and its 
intrinsic bath. If w > wi ow , J\° w (lu) ~ 1/w, correspond- 
ing to 1// noise. 

For a lossy cavity, the bath can be described by a 
Lorentzian spectral density with a central frequency [27| : 



5 2 A 



n[(o 



,) 2 + A 2 



(13) 



which corresponds to a single-mode cavity, with its fre- 
quency broadened by the cavity loss. In Eq. (fT3"|) . A is 
the frequency width of the cavity bath density spectrum, 
w cav is the central frequency of the cavity mode, and g 
denotes the coupling strength between the qubit and the 
cavity. Also, the parameter A is related to the cavity 
bath correlation time and w cav /A is the quality factor Q 
of the cavity. 
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Using Ji{uj) — 9k ~ w k,i), one can derive from 
Eq. © that % is determined self-consistently by the 
equation 



log Vi 



2J i {uj)duj 

(u> + m A) 2 



0. 



(14) 



Below we will solve the equation of motion for the den- 
sity matrix in Hamiltonian (JTJ) and obtain the qubit SE 
spectrum. 



A. Equation of Motion for the Density Matrix 

The equation of motion of the density matrix psb for 
the whole system, i.e., qubit system (S) and bath (B) is 
given by 



Using the Born approximation [28| . the density matrix 
p'g B {t) in Eq. (|2"0j) can be approximated by p'g B {t) = 
p'g(t)p B (0). Tracing over the degrees of freedom of the 
two baths, one obtain 



dt 



p'iit) = -iTr R [V{(t),p s (0)®p B (0)} 



(21) 



-Tr R J [VM, [Vl{t'),p'£{t') ® P b(0)]W, 
a 



where p'g(t) = Ttr [p'sb(^)] ■ Because the two bosonic 
baths are assumed to be in thermal equilibrium, it follows 
from the Bose-Einstein distribution that 



j f PSB(t) = -i[H,p SB (t)}. (15) 
After the unitary transformation (2), we have 

j t p's B (t) = -i[H',p' SB (t)}, (16) 

where p' SB — exp(S)psB exp(— S) is the density matrix 
of the whole system in the Schrodinger picture with the 
transformed Hamiltonian H' [i.e., Eq. ((SJ)]. Below we 
solve the transformed equation of motion, Eq. (|16[) , in the 
interaction picture with Hq — rj Aa z /2+^2 k i u>k,ia\. fli^i- 
In this interaction picture, the transformed Hamiltonian 
H ' can be written as 

V iV) = J2 9k,i a l,i a - ex P [*K,< ~ '? + H - c - ( 17 ) 

The equation of motion for the density matrix p'g B (t) of 
the whole system (S+B) can be written as 



dt 



p'Ut) = -i[vRt), P ^ B (t)} 



(18) 



with p'g B {t) — exp (iHot) p'g B (t) exp (—iHot) . 
Integrating Eq. (|18l) . we have 



Tr 
Tr B 



R (4 
(a k .ia\ 



dk,iPB 
k,iPB 



(22) 
(23) 



where n^^ is the thermal average boson number at mode 
k. Then, substituting Vj{t) into Eq. ([2T]) . we have 



dt 



Ps(t) 



E 

k , i 



~2 

9k,. 



(24) 



f(t') exp [%(u kti - v A)(t - t')\ dt 1 - H.c. 



where 



/ (f) = n k ,i [a^a+p'i (f ) - a + p'i{t')aJ\ (25) 
+ (n fc ,i + 1) \p'i{t')a + a- - a-p'i(t')a+\ . 



PsbW = P'UU) - ^ j M(t%p'L(t')]dt', (19) 

to 

where to is the initial time for the qubit-environment in- 
teraction to turn on. Here we choose to = 0. Substituting 
p's B {t) into Eq. ([18)) , we obtain the equation of motion 
as 



dt 



PsbW - -i[Vi(t), P ' SB (0)] 



(20) 



On the right-hand side of Eq. (|2"4")l . the terms related 
with nk,i and rik,i + 1 describe, respectively, the decay 
and excitation processes, with the rates depending on 
the temperature. Here, for simplicity, we study the zero- 
temperature case with nk.i = 0, i.e., only the spontaneous 
decay occurs, which corresponds to a purely dissipative 
process. The appendix gives the solution of Eq. (|2~4")l 
for the reduced density matrix pg(t) of the qubit in the 
interaction picture. With the solution for p'g(t), one can 
derive the reduced density matrix p' s (t) of the qubit in 
the Schrodinger picture: 
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p' s (t) = exp(ir]Aa z t/2)pg(t)exp(irjAa z t/2) 



(26) 




tAr) £- 



V 



P2l(0) 

P+A+ 



-itArj 



(0) 



p 



Because the reduced density matrix ps(t) in the 
Schrodinger picture with the original Hamiltonian H [i.e. 
Eq. ((TJ] is related to Ps(i) by the relation 

Ps (t) = Tr B [exp (-S) p' s (t)p B exp (S)} . (27) 

Then, using exp (S) = coshl" + c^sinhy, with Y — 

5fc,iCfc,i(fflfc i i-afe,i)/w* i ,i, and tracing over the degrees 
of freedom of the two baths, we obtain 



Psif) 



n 



p' S (t) + ^ L a x p' s (t)a x . (28) 



B. Derivation of the spontaneous Emission 
Spectrum 

When measured by an ideal system with negligible 
bandwidth, the spontaneous emission spectrum can be 
given by 29] 



P{uj) oc Jdt Jdt'exp[-iu;(t-t')]C{t,t'), (29) 



with the two-time correlation function 

C(t,t') = (a+(t)a.(t')) (30) 
= (V<0)k + (*)a_(t')|V>(0)), 

where the c±(t) — exp (iHt) a± exp (—iHt) are the rais- 
ing (+) and lowering (— ) Pauli matrices in the Heisen- 
berg picture and \tp(0)) is the initial state of the whole 
system, which remains unchanged in the Heisenberg pic- 
ture. Using |tA'(0)} = exp (S) \ip(0)) , the two-time corre- 
lation function can be written as 

C(t,t') = (V''(0)|e lff ' t ( 7;e- 4ff '*e lff ' t V^e- iH '*' |^'(0)> , 

(31) 

where a'^_ = exp (S) a± exp (— S). Because the zero- 
temperature case is considered here, this two-time cor- 
relation function can be approximated as 



C(t,t') « (^'(0)|e lff ' t aVe^'*e lff ' t 'a'_e- 4H '*' 



#(0)>, 
(32) 



with 



k 

1+7/ 1—7/ 

= -^r~ a ± H r> a T- 



(33) 



In deriving Eq. (f3"3"|) . the two baths are assumed in the 
ground state for the zero-temperature case. Therefore, 
from Eq. (|33p . we have 



C(t,t') f 



1+7/ 



(a + (t)a- (t')) H , 



«a+ (t)a + {H)) H , + (<r- (t)(7_ (f)>H') 



1-7/ 



(a- (t)o + (t')) H , 



(34) 



where the expectation value of the operator in the trans- 
formed Hamiltonian H ' is 



{(Ta (t) <Tfi (t')) H , 



(35) 



= ^'{Q)\e iH ' t a l a e- iH ' t e iH ' t '(x'pe- iH ' t ' |^'(0)> , 

with a, j3 = ±. 

In our case, 7/ is very close to 1. Therefore, Eq. (|3"4")l 
can be further approximated as 



C{t,t') n C(t,t') H , 



(36) 



with C(t,t')H r = (o+ (t) <t_ (t')) H , . For a qubit state 
specified by a density matrix p(t), we can formulate the 
expectation values of cr+(i), c_(t) and er + (i)<7_(£) by 
the matrix elements, (er + (i)) = (er_(i))* = /?2i(£) and 
C{t,t) = pu(t). According to the quantum regression 
theorem 12911 . the correlation function becomes 



C{t,t + r) H , 



(37) 



-iAr/r I 



In a spontaneous emission process, the initial state is 

an excited state \ip(0)) = exp (-5) |f) ® II l°fc,i> Ofc.a)) 

fe 

which can be achieved by cr^ |<7) . Then, the initial state 
in the transformed Hamiltonian ([5]) is \ip'(0)) = |t) ® 

II |0fc,i,0 fc , 2 ) , i.e. p' n = 1. Therefore, from Eq. (jA~T2|) in 

fe 

the Appendix A, the dynamics evolution of p\ x can be 
expressed as 



P'n(t)=£- 1 



1 



(38) 



Using the Schrodinger equation (see Appendix B) [30| . 
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we have 



1 



p + A+ + A- 



(39) 



P + A 



P + A- 



where £ 1 (^qr^) and £ 1 (^4-) are conjugate 
quantities (see Appendix A). From 

C(t,t + r) H , (40) 

= e -^(r + t) £ -l ( x £-1 ( —L-) e <±nt 
\P+ A +)(r+t) \P+ A -)t 

we obtain the two-time correlation function for any t and 
t' 



C(t,t') H , = £- 



1 



P + A 



-iArjt 



(41) 



+ /t 



-iA-qt' 



Finally, using the Wiener-Khinchin theorem, the SE 
spectrum is given by 



P(w) oc / dt J dt exp \-iu>(t - 1 ) C(t,t) (42) 


1 



T 



£- 



P + A^ 
1 



i — iArjT 



[lj- An- R(lj)} 2 + T{uj) 2 ' 



3, 



c 
o 

— 
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o 
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III. DEPENDENCE OF THE SPONTANEOUS 
EMISSION SPECTRUM ON THE BATHS 

We will show the SE spectrum of the qubit in resonance 
with the cavity central frequency (A = w cav ) as a func- 
tion of the microwave probe frequency for three cases: 
weak, strong and ultra-strong qubit-cavity couplings. 

Weak coupling means the qubit-cavity coupling 
strength g is less than the sum of the dissipation rate 
of the qubit and the cavity. The dissipation rates of the 
qubit due to its intrinsic bath is approximately denoted 
as F q b, which can be A aohm or A ai ow . And the dissi- 
pation rate due to the cavity bath can be approximately 
as the spectrum width of the cavity spectral density A. 
So weak coupling is express as g < (r q b + A) . 

Strong coupling means that the qubit-cavity coupling 
strength g is larger than the sum of the dissipation rate 
of the qubit and the cavity: g > (r q b + A), but it is 
typically two orders of magnitude smaller than the qubit 
energy spacing A and the cavity frequency u; cav , i.e., g ~ 
10 _2 A, such as the case in Ref. fl4l ]. 



FIG. 2: (Color online) The spectrum density of the qubit 
environment, (a) Lorentzian cavity bath and low- frequency 
intrinsic bath of the qubit. (b) Lorentzian cavity bath and 
Ohmic intrinsic bath of the qubit. From Fig. (a) and (b), 
we see that the dominant regimes of the low-frequency and 
Ohmic qubit bath spectral density are different. 



Ultra-strong coupling means that the qubit-cavity cou- 
pling g is to a significant fraction of the transition fre- 
quency A (e.g., g > 0.1A). This case extends to the 
fine-structure limit for the maximal value of an electric- 
dipole coupling. 

The energy spectral density of a bath plays an impor- 
tant role in determining the energy-shift direction and 
the asymmetry of the SE spectrum. In Fig. 2, we show 
the spectral densities for both intrinsic and cavity baths 
of the qubit. For the intrinsic bath of the qubit, both 
a low-frequency bath and an Ohmic bath are consid- 
ered. The spectral density of the cavity bath is symmetric 
about the central frequency aj cav of the cavity. 
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Considering the experimental parameters [lj, |3l| , we 
assume the qubit energy spacing to be A = 10 GHz. The 
dimensionless coupling strength a between the qubit and 
its intrinsic bath (either Ohmic- or low-frequency bath) 
is fixed at a — 10~ 4 , which implies that the decay rate 
of the intrinsic bath is T q b ~ 1 MHz. 



A. Effect of the cavity bath on the spontaneous 
emission spectrum 

To illustrate the effect of the cavity bath with a sym- 
metric spectral density, in Fig. 3 we show the qubit SE 
spectrum when only the cavity bath is present. To en- 
hance these features further, we choose a low quality fac- 
tor Q — 10 2 , and plot the SE spectra in Fig. 3(a) for 
strong (g = 10 2 MHz) and ultrastrong (g = 10 3 , 2 x 10 3 
MHz) qubit-cavity coupling. Due to the scope of ap- 
plication of this method, only the lower bound of the 
ultra-strong coupling g = 0.1 A and g = 0.2A are con- 
sidered. As a comparison, the results under the RWA 
are given in Fig. 3(b). Also, the SE spectra without 
and with RWA for Q = 10 3 are plotted in Fig. 3(c,d). 
In the case of strong qubit-cavity coupling, the two 
peaks of the vacuum Rabi splitting are nearly symmetric 
about uj = A, almost coinciding with the results un- 
der the RWA. When the qubit-cavity coupling increases, 
the height and the position asymmetry (about qubit en- 
ergy spacing A) of the two peaks becomes more appar- 
ent, in sharp contrast to the symmetric SE peaks ob- 
tained under RWA [see Fig. 3(a,b) and (c,d)]. As we 
know, if the RWA is used, the qubit-cavity coupling term 

52Sfe,2( a fe 2 + °fc.2)( cr + + °~-) in the Hamiltonian H be- 

fc 

comes X)#fc-2( a ! 2 a - + a fc,2C+). The energy spectral den- 

k ' ' 

sity in the regions lower and higher than the central fre- 
quency of the cavity (related to absorbing and emitting 
a single photon in the cavity) are identical. Therefore, 
when the qubit energy spacing is in resonance with the 
cavity central frequency, the coupling strength is sym- 
metric about the qubit energy spacing for the absorption 
and emission processes. 

While taking into account the anti-rotating terms, the 

coupling term becomes ^ gk,2{&\ 2 G - + a k,2&+) in the 

k 

transformed Hamiltonian H' , with a renormalized cou- 
pling strength g k>2 = 2i] 2 A# fe>2 / (wjt,2 + n 2 A) . Ob- 
viously, the renormalized coupling strength c/k.2 induces 
the spectral asymmetry: for a symmetric spectral den- 
sity of the cavity bath , in the region ojk.2 < A, due to 
2r]2 A/ (u>k,2 + V2 A) > 1, the renormalized interaction 
9k.2 is larger than gk,2- However, in the region uJk.2 > A, 
owing to 2?7 2 A/ (u;^ 2 + »/2 A) < 1, the effective coupling 
strength gk.2 is smaller than gk,2- 

These results (with and without the RWA) indicate 
that the RWA cannot be used in the range of ultrastrong 
qubit-cavity coupling. The general tendency observed 
here is that the RWA overestimates the frequency shift in 
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FIG. 3: (Color online) Spontaneous emission spectra of the 
qubit only in the cavity bath (symmetric spectral density). 
Figure (a) and (c) are the results without RWA with Q = 
10 2 and Q = 10 3 . Fi gures (b) and (d) are the results of 
the RWA with Q = 10 2 and Q = 10 3 . To see the height 
asymmetry of two peaks clearly, the horizontal grid lines are 
plotted as reference. Note that in (a) and (c), the two peaks 
of spontaneous emission spectrum present obvious height and 
position asymmetry (about ui — A) in ultra-strong qubit- 
cavity coupling. 



the low-energy regime uj ~ —g, while it under- estimates 
the frequency shift in the high-energy regime ui ~ g. Our 
results are consistent with the results in Ref. flTTl . 



B. Combined effect of both intrinsic and cavity 
baths on the spontaneous emission spectrum 

Although a high-Q seems plausible for minimizing the 
loss of the cavity, it limits the measurement speed. Here, 
we consider a cavity with the quality factor [9|, [l4| Q = 
10 4 in the presence of a qubit (see Fig. 4). We also plot 
the SE spectrum for Q = 10 3 in Fig. 5, and see how the 
quality factor affects the results. The dissipation rate of 
the cavity is approximately to the spectral width A of the 
cavity bath. If Q = 10 4 , the dissipation rate of the cavity 
bath is about 1 MHz (the same order of magnitude of the 
bath dissipation rate). Figures 4(a) and 4(b) show the 
spectra of the qubit coupled with a low-frequency and 
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FIG. 4: (Color online) Spontaneous emission spectra of the 
qubit in resonance with the central frequency o; cav of the 
cavity for weak, strong and the lower bound of the ultra- 
strong qubit-cavity interactions. The cavity quality factor 
is Q — 10 4 . (a) coupling strength to the low-frequency in- 
trinsic bath of the qubit ai ow = 10 -4 , (b) coupling strength 
to the Ohmic intrinsic bath of the qubit aohm = 10~ 4 . To 
see the height asymmetry of the peaks clearly, the horizon- 
tal grid lines are plotted as reference. Note that (a) demon- 
strates the obvious height asymmetry in the case of strong 
and ultra-strong qubit-cavity coupling and the asymmetry in- 
creases as the qubit-cavity coupling grows. Figure (b) shows 
inverted height asymmetry of two peaks from (a) in the strong 
qubit-cavity coupling case, but as the qubit-cavity coupling 
increases to the ultra-strong regime, the height asymmetry of 
the right and left spectral peaks are inverted. 
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FIG. 5: (Color online) Spontaneous emission spectra of the 
qubit in resonance with the central frequency oj cav of cavity for 
weak, strong and the lower bound of the ultra-strong qubit- 
cavity interactions. The cavity quality factor is Q = 10 3 . (a) 
coupling strength to the low-frequency intrinsic bath of the 
qubit ai ow = 10 -4 , (b) coupling strength to the Ohmic intrin- 
sic bath of the qubit aohm = 10 -4 . To see clearly the height 
asymmetry of the peaks, horizontal grid lines are plotted 
as reference. Note that (a) demonstrates an obvious height 
asymmetry in the case of strong and ultra-strong qubit-cavity 
coupling and the asymmetry increases as the qubit-cavity cou- 
pling grows. Figure (b) shows the inverted height asymmetry 
of two peaks from (a) in the strong qubit-cavity coupling case, 
but as the qubit-cavity coupling increases to the ultra-strong 
regime, the height asymmetry of the right and left spectral 
peaks are inverted. 



an Ohmic intrinsic baths, respectively. From Fig. 4(a), 
we see that in the case of weak qubit-cavity coupling 
g = 10~ 5 A = 0.1 MHz, the SE spectrum is a single peak 
with the central frequency larger than the energy spacing 
A of the bare qubit, which corresponds to a blue-shift. In 
the case of strong qubit-cavity coupling, g = 2 x 10 _2 A = 
2 x 10 2 MHz, the SE spectrum shows the vacuum Rabi 
splitting, with the two height asymmetric peaks, just as 



shown in Rcf . Il4|. By further increasing the qubit-cavity 
coupling g = 10 _1 A = 1 GHz, not only the height of the 
SE peaks but also their positions demonstrate a strong 
asymmetry about w = A. The qubit-cavity coupling 
strength g = 10 _1 A is the lower bound of the ultra- 
strong coupling regime. 

Figure 4(b) shows the SE spectrum of a qubit in the 
Ohmic intrinsic bath. For a weak qubit-cavity coupling, 
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the central frequency of the SE spectrum shifts to an en- 
ergy slightly lower than the energy spacing A of the bare 
qubit, which corresponds to a red-shift. This energy-shift 
direction is opposite to the case when the qubit is in a 
low- frequency intrinsic bath [see Fig. 4(a)], and the en- 
ergy shift is smaller. For a strong qubit-cavity coupling, 
the two SE peaks also show a weak asymmetry , with the 
left peak higher than the right peak. This peak asym- 
metry is inverted from the same case in Fig. 4(a). As 
the qubit-cavity coupling increases to the ultra-strong 
regime, the height asymmetry of the left and right SE 
spectrum are inverted [see Fig. 4(b)]. The SE spectra 
of the qubit in a cavity with quality factor Q = 10 3 are 
shown in Fig. 5, which present nearly the same features as 
in Fig. 4, in addition to the broader SE peaks in Fig. 5. 
This is because of an increased dissipation rate of the 
qubit induced by a larger cavity dissipation rate. These 
results are briefly summarized in table I. 

Figures 4 and 5 show that the two SE peaks of the 
qubit in both a low-frequency and Ohmic intrinsic baths 
have very different behaviors; the right SE peak is higher 
than the left SE peak in both strong and ultrastrong 
qubit-cavity coupling regimes, while the right SE peak 
is lower than the left SE peak in the strong qubit-cavity 
coupling regime, and in the lower bound of the ultra- 
strong qubit-cavity coupling regime, the left SE peak is 
slightly lower than the right peak. Therefore, the differ- 
ent asymmetric behaviors of the SE spectrum of the qubit 
in the low-frequency and Ohmic baths may be used to 
distinguish the intrinsic noise of the qubit. In the exper- 
iment in Ref. [T3 |. the height asymmetric SE spectrum 
of the qubit in the strong qubit-cavity coupling regime 
shows that the right SE peak is higher than the left SE 
peak. This reveals that the low-frequency intrinsic noise 
is dominant in the superconducting qubit in Ref. \lA} . 



IV. DISCUSSION AND CONCLUSION 

Below we discuss the reason why the SE spectrum in 
the low-frequency intrinsic bath is different from that in 
the Ohmic intrinsic bath. We begin with the standard 
Born-Markov master equation (at zero temperature) 0, 
|32| , which is usually derived under the RWA [33[ 

P = -i [H,p] +7(2(T_p(T + - o- + (7_p- pcr+cr-) (43) 
+n (2apa) — a^ap — pa'aj . 

Equation (|43|) is typically used to describe the whole sys- 
tem consisting of an atom and a (single mode) cavity 
with dissipation via two channels: the dissipation of the 
qubit due to a free-space mode (the term proportional to 
7), and the dissipation of the qubit due to the cavity loss 
(the term proportional to k). In quantum optics, the SE 
spectrum is often calculated by substituting the Jaynes- 
Cummings model Hamiltonian H^ c = g (er_a^ + cr+ft) 
as H, into Eq. (|43l) . Then the vacuum Rabi frequency 
spitting with two symmetric SE peaks is derived, in the 



strong qubit-cavity coupling, where the position of the 
two peaks are exactly at to = ±g. This RWA result 
is different from the experimental result |l4T |. where the 
two SE peaks are asymmetric about to — A (the energy 
level spacing of the bare qubit). As shown in Sec. IIIIB1 
the height asymmetry of the SE spectrum is obtained 
beyond RWA. Therefore, it can be concluded that the 
anti-rotating terms produce an asymmetric SE spectrum 
of the qubit. 

In conclusion, we discussed the SE spectrum of a qubit 
in the environment described by two baths: intrinsic bath 
and cavity bath. We only consider that the central fre- 
quency of the cavity mode is resonant with the qubit 
energy spacing (w cav = A). We analyze in detail the 
qubit 's SE spectrum in the weak, strong, and the lower 
bound of the ultra-strong coupling regimes, and compare 
the SE spectra in two kinds of qubit baths: low-frequency 
and Ohmic baths. In the low-frequency bath, the height 
asymmetry of the vacuum Rabi splitting peaks increases 
as the coupling strength grows. However, for the Ohmic 
bath, the height asymmetry of the SE spectrum is re- 
duced, then the height asymmetry of the left and right 
peaks are inverted, when the coupling strength is in- 
creased. All of these results show that for strong qubit- 
cavity coupling, the asymmetry of the splitting peaks in 
the SE spectrum comes from the anti-rotating terms and 
the non-constant spectral density of the bath. 
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Appendix A: Solution of the equation of motion of 
the density matrix 

This appendix offers detailed calculations for solving 
the master equation in Eq. (|2"4"1) . We use the basis |1) = 
\i) and |2> = |t), where a z \i) = ||) , a, |t) = - |t) , 
to define the reduced density matrices of the qubit. By 
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TABLE I: Summary of our main results for the spontaneous emission spectra in the case of weak, strong and ultra-strong qubit- 
cavity coupling. The spectra with symmetric peaks (S), height asymmetric peaks (AS) and very asymmetric peaks (VAS) are 
abbreviated as S, AS and VAS. The asterisks indicate which are inverted peak height asymmetry from the AS in the same case 
of the low-frequency bath. These results are described in detail in the main text. 



Bath 


Cavity 
quality factor 


qubit-cavity coupling 


weak 


strong 


ultra-strong 


Low-frequency 


high Q 


single 
peak 


AS 


VAS 


low Q 


single 
peak 


AS 


VAS 


Ohmic 


high Q 


single 
peak 


AS* 


AS 


low Q 


single 
peak 


AS* 


AS 



using the Laplace transform 

oc 

p(p) = C [p(t)\ = J dtp(t) exp (-pi) 
o 

and the convolution theorem 
t 

C I dt'h{t')f 2 {t-t') 



(A.l) 



C[h(t)]C[h(t)], (A.2) 



the master equation (1241) for the qubit system can be 
solved as 

PPs(p)-Ps(Q) (A.3) 
= -A-a + a-p s (p) + (A+ + A-)a-p s {p)a+ 



where 



(A.4) 



This equation is a Lyapunov matrix equation. 

The Kronecker product property in matrix theory 
shows that 

Vec (M 1 pM 2 ) = Mi ® Ml Vec (p) , (A.5) 

where Vec (p) represents the vector expanding of matrix 
p along rows, and the superscript T denotes the trans- 
pose of the matrix. We expand the matrix equation in 
Eq. (|A.3|) into vectors along rows: 



U (p) Vec p s (p) 



Vec 



Ps(0) 



where 



U(p) = ph + A+I 2 <8> {(T+v-f 
-(A- + A + )o- ®cr£ 

+A_(cr+(7_) ® J 2) 



(A.6) 
(A.7) 



where I n is the n x n identity matrix. Thus, the 2x2 
matrix equation (|A.3[) is transformed to the 4x4 vec- 
tor equation in (IA.6I) . The solution of Eq. (|A.6|) can be 
formally written as 



Vec 



P S (P)] ^Uip)- 1 Vec \p' s (0) 



(A. 



with 



u&r 1 = 



P+A++A. 









P+A 





\ 





p+A + +A 

By using the inverse Laplace transform 



P+A+ 

i / 

P / 



(A.9) 



t(t) = C- 1 [a(p)} = -L f dpa(p)eiq>(pt), (A.10) 
2m J 



we obtain, 

Vec^' (t)] - C^Uip)- 1 Vec [ Ps (0) 
Then, p s (t) is given by 
Ps(t) 



(A.11) 



(A.12) 



Pii(O) 



P+A++A- 

P2l(0) 

P+A- 



Pii(O) 

p 



£12(0) 
P+A+ 

Pii(O) 



P+A++A- 



P22(0) 

P 



Below we calculate the inverse Laplace transform 

^ X (i^l) and ^ (pTaT+a-)- From E q- 
we have 



1 



P + A_ 



1 
2~7ri 



exp (pi) 



p + E^Li/b-iK^-^A)] 

fc,z 

(A.13) 



dp. 
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With p replaced by iuj + + [3J|, the above expression where p stands for the Cauchy principal value. Let 

R(u) = + R 2 (uj) and T{uj) = T^w) + V 2 (u)), then 

we have 



becomes 

c- 1 



1 



P + A- 



exp(iwt) 



(A.14) 



-duo. 



2ni J jj - Y^9k,il I u ' _ ~ ■ ~ ,{)J[ 

k,i 



For the term X)fefffei/( w — u k,i — *0 + ), we denote the 
real and imaginary parts as Ri(ui) and Ti(u)), where i = 
1 for the intrinsic bath and i = 2 for the cavity bath. 
Explicitly, we can write 



1 



1 

2nd 



p + A- 

+ 00 



(A.17) 



exp(iwt) 



lu - R{lo + t] A) - iT(w + r/A) 



du. 



RM = p £-z 



: ,1 



and 



(A. 15) Similarly, £ and C 1 p+A ^ +A _ can also be de- 

rived as 



2t7i A 



Ji(w 



u/ + iji A / (u; — w') 



(A.16) 



1 



P + A- 



(A.18) 



2??, A 
cj + rji A 



and 



P + A^ 



A. 



(A.19) 



exp(iwi) 



2ni J u- R(lo + At?) + i?(Ar7 - u) - i [T(lo + At?) + T(A?/ - w)] 



Appendix B: Solution of the Schrodinger equation 

Below, we will solve the equation of motion of 
wave function beyond the RWA in the transformed 
Hamiltonian W in Eq. ([5]). Since the total exci- 
tation number operator of the qubit-cavity system, 

N = J2k ( a ! a fc + bi-bkj + (1 + &z) /2 in the transformed 
Hamiltonian is a conserved observable, i.e., [N,H'\ = 0, 
it is reasonable to restrict our discussion in the single- 
particle excitation subspace. A general state in this sub- 
space can be written as 

i*(*)> = x (t) \s2) n \o k ,iOk,2)+Y,pkM i*i> n I *,! 1 *.* 



we have 



k,i 



(B.l) 

k jl k ,i}) means either cavity bath or 



where the state 

qubit spontaneous dissipation bath with one quantum ex- 
citation. Substituting |$(t)) into Schrodinger equation, 



••^ = ^x(*)+E^ w (t) 1 

k.i 



dt 



I = (Wfc,- 



^)^,i(*)+E^x(*). 



dt v ' 2 
Applying the transformation 

X {t) = X{t) exp ( -i^t 



0k,i(t) = /3 fc;i (t)exp 



-H W M Y' 



Eqs. (IB. 21) and (|B.3|) is simplified as 

= — '£ 



(B.2) 
(B.3) 

(B.4) 
(B.5) 



= -i E Vk,iPu,i (*) exp [-i(wfe,i - ?7 > ( B -' 
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dt 



= -iVfe,iX(*) exp [i(u kyi - 77 A)t] . (B.7) 



Integrating Eq. (|B.7[) and substituting it into Eq. (|B.6[) . 
we obtain 



= -/ E v M ex PHKi-'? A )(*-* / )]x(*')*'- 
*m 

(B.8) 



This integro-differential equation (|B.8|) is solved exactly 
by Laplace transformation, 



x(p) 



x(o) 



fc.t 



(B.9) 



When the initial state is an excited state \ip'(0)) — 
It) ® II |0fe,i,0 fe , 2 ) , i.e., x(0) = 1. Applying the Inverse 



Laplace transformation, we get 



m = c 



P + A 



+ /t 



(B.10) 



From Eq. (|A.18|) . the dynamics evolution of p' n can be 
expressed as 



Pn(t) = X*(t)x x (t) 



(B.ll) 



P + A 



x C- 1 



P + A- 
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